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Abstract 

A flight test was performed with the NASA 
Lewis Research Center's DH-6 icing research air- 
craft. The purpose was to employ a flight test 
c3 procedure and data analysis method, to determine 
the accuracy with which the effects of ice on air- 
uj craft stability and control could be measured. 

For simplicity, flight testing was restricted to 
the short period longitudinal mode. 

Two flights were flown in a clean (baseline) 
configuration, and two flights were flown with 
simulated horizontal tail ice. Forty-five repeat 
doublet maneuvers were performed in each of four 
test configurations, at a given trim speed, to 
determine the ensemble variation of the estimated 
stability and control derivatives. Additional 
maneuvers were also performed in each configura- 
tion, to determine the variation in the longitudi- 
nal derivative estimates over a wide range of trim 
speeds . 

Stability and control derivatives were esti- 
mated by a Modified Stepwise Regression (MSR) 
technique. A measure of the confidence in the 
derivative estimates was obtained by comparing the 
standard error for the ensemble of repeat maneu- 
vers, to the average of the estimated standard 
errors predicted by the MSR program. A multiplic- 
ative relationship was determined between the 
ensemble standard error, and the averaged program 
standard errors. In addition, a 95 percent confi- 
dence interval analysis was performed for the ele- 
vator effectiveness estimates, C m&e . This analysis 
identified the speed range where changes in C mse 
could be attributed to icing effects. 

The magnitude of icing effects on the deriva- 
tive estimates were strongly dependent on flight 
speed and aircraft wing flap configurai ton. With 


*3oint Institute for Advancement of Flight 
Science, NASA Langley Research Center. 


wing flaps up, the estimated derivatives were 
degraded most at lower speeds corresponding to 
that configuration. With wing flaps extended to 
10°, the estimated derivatives were degraded most 
at the higher corresponding speeds. 

The effects of icing on the changes in longi- 
tudinal stability and control derivatives were 
adequately determined by the flight test procedure 
and the MSR analysis method discussed herein. 

List of Symbols and Abbreviations 
b wing span, m 

C],C m , and C n roll, pitch, and yaw moment coefi- 
cients, respectively 


Cy ,Cy , and Ci longitudinal, lateral, and verti- 
cal force coefficients, 
respectively 


c 

mean aerodynamic chord, m 

c.g. 

airplane center of gravity 

g 

acceleration due to gravity 

Ix,Iy, and I i 

moments of inertia, kg-m 2 

l XZ 

product of inertia, kg-m 2 

JSR 

cost function for Stepwise Regres- 
sion Algorithm 

KIAS 

knots indicated airspeed, kts 

m 

aircraft mass, kg 

N , n 

number of points 


roll, pitch, and yaw rates, respec- 
tively, rad/sec or deg/sec 


This paper is declared a work of the U.S. Government and is 
not subject to copyright protection in the United States. 
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wing area, m 2 
shaft horsepower, hp 


small ice accretions and to those larger accre- 
tions caused by ice protection systems failures 
or operator errors . 


longitudinal, lateral, and verti- 
cal velocity components 

total airspeed, m/s or knots as 
given by (u 2 + v 2 + w 2 ) 1 ^ 2 

measured aircraft response or con- 
trol surface input 

linear approximation of C m&e corre- 
sponding to some x G (velocity) 

computed aerodynamic coefficients 

corrected (free stream) angle of 
attack, rad or deg 

angle of sideslip, rad or deg 

elevator deflection, rad or deg 

aerodynamic stability or control 
parameters 

constant offset term 

derivative w.r.t. some response or 
control surface input 

pitch and roll angles, respec- 
tively, rad or deg 

air density, kg/m 2 


variance estimate of measurement 
noi se 


Introduction 


Currently, no analysis methodologies exist 
that provide the capability to predict the effect 
of icing on aircraft flight characteristics. 
Therefore, within the scope of the NASA Lewis air- 
craft icing research program, a high priority has 
been assigned in measuring the effects of icing on 
aircraft performance and stability and control. 
Experimental data acquired from flight tests pro- 
vide the basis for validating icing effects meth- 
odologies, and for establishing correlations 
between computational codes, wind tunnel tests, 
and full scale flight tests. Performance test- 
ing, which was extensively discussed in Refs. 1 
and 2, provided reasonably consistent results and 
a high degree of confidence in measurement accu- 
racy. As a result, these data are now being used 
in validating current methodologies that predict 
ice shape growth and their effect on aerodynamic 
performance. However, initial flight testing to 
estimate effects of icing on aircraft stability 
and control derivatives, showed that the results 
could be obscured by the inherent problem of not 
knowing the true accuracy of the derivative esti- 
mates. Stability and control flight testing as 
reported in Refs. 3 and 4 used an early version 
of a Modified Maximum Likelihood Estimation (MMLE) 
technique to estimate the derivative values. The 
accuracies of these derivative estimates were 
obtained by an approximation to the Cramer-Rao 
bound given by the maximum likelihood algorithm. 
Because it was known that these error estimates 
indicated higher accuracies than those which can 
be achieved by repeated experiments, 5 and because 
the differences between derivatives from the base- 
line and iced aircraft were found to be small, it 
was decided to conduct further flight testing to 
determine a realistic accuracy with which stability 
and control derivatives could be estimated. The 
purpose of this paper is to report on the results 
of those tests, and to present the conclusions 
drawn from the data analysis programs. 


Modern transport aircraft, certified for 
flight in icing conditions, normally employ ice 
protection systems to prevent ice from accumulat- 
ing on critical lifting components. If these 
systems fail, or are improperly operated, ice can 
form on these components and pose a serious threat 
to safe flight characteristics. This problem, 
which is a continual concern for current genera- 
tion aircraft, will be of even greater concern for 
future generation aircraft employing advanced pro- 
pulsion systems and optimized aerodynamic designs. 
For example, high propulsive efficiency engines 
such as propfans, may not provide sufficient bleed 
air for wing and tail anti-icing. The use of 
energy-efficient de-icing systems, which require 
some small build-up of ice before activation, will 
be their likely replacement. Aircraft equipped 
with these systems will have to demonstrate 
acceptable flying qualities with some leading edge 
ice contamination, especially when flown in the 
take-off or landing configuration. This require- 
ment also applies to aircraft whose aerodynamic 
efficiency, enhanced by advanced airfoil designs 
and relaxed static stability criterion, is degraded 
when small leading edge ice formations occur 
between de-icing cycles. Likewise, aircraft 
employing full flight control augmentation must 
also display an acceptable tolerance to these 


For simplicity, the flight test plan was 
structured to examine only the longitudinal short 
period characteristics of the research aircraft. 
Four data collection flights were defined: two 

flights, were flown clean; and two flights were 
flown with an artificial ice shape (described later 
in this report), attached to the leading edge of 
the horizontal tai 1 . 

Flight test maneuvers consisted of numerous 
repeat elevator doublets at a specified trim con- 
dition to determine the ensemble standard devia- 
tion for the parameter estimates. Additional 
maneuvers were also flown over the entire level 
flight angle of attack envelope to obtain deriva- 
tive estimates for different flight conditions. 
Flight test maneuvers were analyzed using a modi- 
fied stepwide regression (MSR) algorithm. The 
results are discussed relative to the accuracy of 
the derivative estimates, and to the differences 
between the estimates obtained for the uniced 
versus iced horizontal tail. Some conclusions 
are also drawn as to the physical interpretation 
of the derivative estimates and their relevance 
to aircraft handling characteristics. Following 
the introduction, this paper is organized in sec- 
tions which describe the research airplane. 
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instrumentation, flight test maneuvers, data 
reduction and analysis techniques, accuracy of 
the derivative estimates, effects of tail ice on 
handling characteristics and conclusions. 

Research Aircraft 

The NASA Lewis icing research aircraft shown 
in Fig. 1, is a modified DeHavilland DH-6 Twin 
Otter. It is powered by two 550 SHP Pratt and 
Whitney PT6A-20A turbine engines driving three 
bladed Hartzell constant speed propellers. Physi- 
cal dimensions, mass, and inertial specifications 
for the aircraft are found in Table I. 

The flight controls are mechanically operated 
through a system of cables and pulleys. Ailerons 
droop as wing flaps extend. Their degree of move- 
ment increases in proportion to the amount of wing 
flap extension. There are two trim tabs on the 
elevator. The right tab is interconnected with the 
wing flap system to provide automatic longitudinal 
compensation as the wing flaps are operated. The 
left trim tab is an adjustable tab, which is oper- 
ated by the pilot. 

For flights where an artificial ice shape was 
used to simulate a moderate glaze icing condition 
on the tail, formed aluminum caps were placed over 
the left and right leading edges of the horizontal 
tail, and held in place with straps and clamps. 
Artificial ice cut from styrofoam was then secured 
to the tail caps with double sided tape. The 
geometry and positioning of the artificial ice was 
based on an actual tail icing condition recorded 
on an earlier research flight in natural icing 
clouds. The entire arrangement is shown in Fig. 2. 

Instrumentation 

The icing research aircraft was equipped with 
three linear control position transducers (CPT) to 
measure the position of the elevator, ailerons, 
and rudder. The elevator CPT measured elevator 
position at the control horn to eliminate correc- 
tions for cable stretch. Flap position, constant 
for a given maneuver, was manually recorded from 
the flap position indicator in the cockpit. 

Airspeed, pressure altitude, angle of attack 
and angle of sideslip were measured by means of 
sensors located on the flight test nose boom, 

Fig. 3. Pitot-static pressures were sensed by a 
Rosemount 858 probe, which had been calibrated by 
the "trailing cone" method. The high frequency 
response angle of attack and angle of sideslip 
vanes, were constructed from balsa wood, and bal- 
anced by a "heavy-metal" counterweight. Angle of 
attack was referenced to the fuselage reference 
line, which is coincident with the aircraft floor 
line, and angle of sideslip was referenced to the 
aircraft longitudinal axis. 

Kohlman Systems Research, Inc., Lawrence 
Kansas under contract to NASA, provided a stability 
and control data acquisition system (DAS). This 
system consisted of an attitude gyro, three axis 
orthogonal accelerometers , three axis angular rate 
gyros, and a digital data recording system which 
recorded 21 channels of flight test parameters 
with 12 bits per channel. All data were recorded 
on-board at 20 Samples per Second (SPS), and all 
channels were sampled within one millisecond by 


the DAS. Because of buffer limitations, only 
24 sec of 20 SPS data could be obtained before the 
buffer contents had to be written on to a cas- 
sette tape. The write-off period was nominally 60 
to 90 sec. Hence, all flight test maneuvers were 
planned to last no more than 24 sec. During the 
write-off period, the aircraft would be reposi- 
tioned and trimmed for the next data point. The 
list of recorded parameters, ranges, and sensitiv- 
ities is provided in Table II. 

FI iqht Testing 

The center of gravity of the aircraft was 
measured experimentally. These measurements 
represented the zero-fuel mass characteristics of 
the aircraft and were corrected in the analysis 
program for actual crew and fuel weights recorded 
during each flight test maneuver. 

Longitudinal stability and control flight 
testing was accomplished in four airplane configu- 
rations defined according to wing flap setting and 
the simulated icing condition of the horizontal 
tai 1 plane. They were : 

(1) Clean tail, wing flaps up (0°) 

(2) Clean tail , wing flaps 10° 

(3) Artificial horizontal tail icing, wing 
flaps up (0°) 

(4) Artificial horizontal tail icing, wing 
flaps 10° 

A constant engine power setting associated 
with each wing flap condition, was used while 
executing flight test maneuvers. This eliminated 
engine power as a variable in the derivative esti- 
mates. Two hundred seventy-five SHP/engine was 
used for the 0° wing flap cases, and 219 SHP/ 
engine was used in the 10° wing flap cases. 

The standard small amplitude input was a 
"2-1" elevator doublet which was designed to 
excite the short period longitudinal mode of the 
airplane. Elevator control column deflection was 
input by the pilot to effect a planned ±0.3 g ver- 
tical acceleration from 1.0 g level flight, how- 
ever, due to an erroneous cockpit accelerometer, 
these maneuvers were more on the order of ±1 .0 g 
vertical acceleration. A typical example of one 
of these doublet maneuvers is shown in Fig. 4. 

Forty-five sets of identical small amplitude 
doublets were performed at a given flight condi- 
tion in each configuration to determine the accu- 
racy of the longitudinal stability and control 
derivative estimates. For the 0° flap configura- 
tions, repeat doublets were initiated from a 
level, unaccelerated flight condition at a trim 
speed of 120 KIAS. For the 10° flap configura- 
tions, these doublets were initiated from a 
level, unaccelerated flight condition at a trim 
speed of 100 KIAS. 

Additional small amplitude doublets were 
then performed in each configuration at incremen- 
tally higher angles of attack (lower trim speeds). 
Power settings associated with each flap configu- 
ration remained constant for these maneuvers, 
therefore doublet maneuvers initiated at the lower 
trim speeds were performed in climbs through a 
median test altitude. 
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Data Reduction and Analysis Technique 

The recorded casette tapes from each flight 
were initially sent to KSR, Inc., where the data 
were reduced to engineering units and rewritten 
on nine-track magnetic tape. These data sets were 
then transferred to the NASA Langley Research Cen- 
ter, Aircraft Guidance and Control Branch, for 
stability and control analysis. 

The first step in the analysis procedure was 
to correct the measured data for all known errors. 
All sensor data were corrected for c.g. offset, 
angle of attack data were corrected for upwash and 
boom bending effects, and several flight maneuvers 
were analyzed with a maximum liklehood algorithm 6 
to estimate instrument biases and scale factors. 
These estimated corrections, when applied to the 
measured data, produced data sets which were ready 
for model structure determination and parameter 
estimation. 

For the model structure determination proce- 
dure, it was assumed that, 

(1) The general equations of rigid body 
motion adequately define the airplane motion. 

(See appendix) 

(2) The model for the aerodynamic force-and- 
moment-coeff icients can be represented by multi- 
variable polynomials in response and control 
variables. The parameters in these equations are 
the coefficients of the Taylor series expansion 
around the values corresponding to the initial 
steady-state flight. 

(3) Linear terms in the Taylor series expan- 
sion make generally larger contribution to aerody- 
namic functions, followed by higher-order terms. 

The third assumption, resulting in a con- 
straint on the selection of significant terms in 
the regression equation, is what gives rise to the 
Modified Stepwise Regression (MSR) technique. 

MSR is a modified version of linear regres- 
sion, which can determine the structure of the 
aerodynamic model equations, and estimate the 
model parameters. The determination of an ade- 
quate model for the aerodynamic coefficients 
includes three steps: the postulation of terms 

that might enter the model, selection of an ade- 
quate model, and the verification of the model 
selected. The general form of aerodynamic model 
equations can be written as 

y( t) - e 0 + Oi X] (t) + . . . + 0 n x n (t> (1) 

where y(t) represents the resultant coefficient 
of aerodynamic force or moment. In the polynomi- 
nal representation of the aerodyamic coefficient, 
01 to 0 n are the stability and control deriva- 
tives, 0 O is the value of any particular coeffi- 
cient corresponding to the initial steady-flight 
conditions, and xi to x n are the regressors 
formed by the airplane output and control varia- 
bles or their combinations. 

For example, if we let y - C m , the airplane 
pitching moment coefficient, we may choose the 
candidate variables, X], X£,. . . X n , as: 
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and their corresponding coefficients as the 
stability and control derivatives: 
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After postulating the aerodynamic model equa- 
tions, the determination of significant terms 
among the candidate variables and estimation of 
the corresponding parameters follows. The varia- 
ble chosen for entry into the regression equation 
is the one that has the largest correlation with 
y after adjusting for the effect on y of the 
variables already selected. The parameters are 
estimated by minimizing the cost function 


J 


SR ~ 


N 


£ 

1-1 



a ~ 12 

-E 


M 


(4) 


where N is the number of data points and (ft, + 1) 
is the number of parameters in the regression 
equation at any step. 


At every step of the regression, the varia- 
bles incorporated into the model in previous stages 
and a new variable entering the model are reex- 
amined. Any variable that provides a nonsignifi- 
cant contribution (due to correlation with more 
recently added terms) is removed from the model. 

The process of selecting and checking variables 
continues until no more variables are admitted to 
the equation and no more are rejected. Experience 
shows, however, that the model based only on the 
significance of individual parameters in the model 
in Eq. (1) can still include too many terms and, 
therefore, may have poor prediction capabilities. 
Several criteria for the selection of an adequate 
model are introduced in Ref. 7, and the details of 
the entire procedure are explained in that 
reference. 


Applying these criteria, an adequate model 
equation for the airplane pitching moment coeffi- 
cient was defined by a Taylor Series expansion, 



a + 
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( 5 ) 


the solution of which provided an adequate fit to 
the measured flight data. In a similar manner, 
this procedure was used to estimate the other lon- 
gitudinal stability and control derivatives 
described in this report. The analyses for the 
C z and C m derivative estimates, and the rela- 
tive effects of tail icing on these estimates is 
given in Figs. 5 to 8. A physical interpretation 
of these results will be discussed in a later 
section of this report. 
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Accuracy of the Derivative Estimates 
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The 45 repeat doublet maneuvers, identically 
executed in each of the four configurations, pro- 
vided a statistically significant data base for 
parameter and error estimate accuracy studies. 

The standard error calculated for each ensemble 
of parameter estimates verified the effectiveness 
of maneuver selection through demonstrated repeat- 
ability of the predicted parameters. In addition, 
a comparison of the average of the standard error 
values estimated by the computer program to the 
standard error for the ensemble, provided an eval- 
uation of how realistic the program estimated 
standard error values were. 

Repeat maneuvers were analyzed using the Mod- 
ified Stepwise Regression and the resulting mean 
values of coefficients were compiled in Table III. 
Accompanying each parameter estimate is the aver- 
age standard error associated with it as given by 
the estimation algorithm, the standard error for 
the ensemble of values from the 45 repeat maneu- 
vers, and the ratio of the two. The ratio given 
for each parameter is the factor by which the pro- 
gram estimated standard error for that parameter 
should be multiplied to produce a realistic 
expected standard error. For example, the stand- 
ard error estimated by MSR for C z should be 
multiplied by 2.9 to achieve a reaTistic approxi- 
mation to an actual ensemble standard error. 


A fundamental issue throughout this research 
is that of distinguishing between the effects of 
ice and the inherent variability of predicted 
parameters. In addition to the error analysis 
performed on the values obtained from the 45 
repeat maneuvers, the C m&e estimates for the 0° 
flap condition were considered (Fig. 6(c)). This 
set of points was chosen for this portion of the 
analysis because each of the two groups of values 
(iced and uniced) lend themselves to linear 
approximations. Also, determining a prediction 
interval in this particular case serves to quan- 
titatively enforce what is qualitatively a very 
noticeable difference in parameters for the iced/ 
uniced cases. After approximating the two groups 
of points by straight lines, prediction intervals 
about these lines were calculated as follows: 


W + ^*'2, n-2 


where 


+ n 


-1 


+ Sxx J 


( 6 ) 


y 0 (x 0 ) linear approximation of 

C m corresponding to 
some (velocity) value 

a* 0.05 to achieve a 95 percent 

prediction interval 



standard error for values of 



s„(y) * Sp A <y) ensemble to program 

average obtained from 
Table III. 


S * (x. - x; sum of squares of residuals 

XX 1 

The resulting interval is one within which we 
can expect 95 percent of all predicted elevator 
effectiveness values to fall. Figure 9 presents 
these prediction intervals and gives an illustra- 
tion of both the areas in which icing effects are 
discernable as well as the speed range in which we 
cannot attribute any change in C m&? to the arti- 
ficial ice shape. A similar analysis can be done 
on other parameters by fitting them in a piecewise 
linear fashion. 

Effects of Tail Ice on Handling Characteristics 

The effects of icing on the horizontal tail 
plane were obtained by analyzing all of the dou- 
blet maneuvers executed throughout the attainable 
ranges of speed and angle of attack. The MSR ana- 
lyzed results of these maneuvers, in the form of 
the estimated stability and control derivatives, 
are presented in Figs. 5 to 8. Each parameter 
value is plotted with ±2a error bars based on the 
standard error as given by the computer program. 

The MSR analysis indicated that no nonlinear terms 
were required for an adequate aerodynamic model. 

Figures 5 and 6 show an analysis for the 0° 
flap cases. The C z derivatives indicate an 
apparent loss of elevator effectiveness, especially 
at low airspeeds that correspond with higher 
angles of attack. The pitching moment derivatives, 
Cm a » Qnq* anc * C m g e , which have an important 
effect on airplane short period stability and con- 
trol characteristics, are also more strongly 
affected by the ice contaminated horizontal tail 
at low speeds. At higher speeds, the effects of 
tail ice become negligable. This analysis was 
supported by pilot comments reporting that at low 
speeds, the "iced" airplane was less responsive 
to elevator control inputs, and more weakly damped 
in pitch. These comments, however, were not based 
on specific piloting tasks normally used in defin- 
ing pilot ratings, but on the perceived relative 
differences (clean versus iced) in longitudinal 
control force, response, and pitching moment cha- 
racteristics while performing parameter identifi- 
cation maneuvers. Further, the magnitude of these 
effects were not judged to be great enough to 
change pilot ratings if such an evaluation were 
to have been carried out. 


n number of points 

*<x*/2,n-2 student t-di str i bution 

s standard error (can be one 

of three listed below) 

sp^(y) * average of standard errors obtained 
from the computer program 


In the ten degree flap cases, Figs. 7 and 8, 
the effects of horizontal tail icing show an oppo- 
site trend when compared with the 0° cases. 
Analysis of the stability and control derivatives 
show slightly degraded pitching moment coeffici- 
ents up to speeds of approximately 90 kn. Above 
these speeds, there is a more pronounced degrada- 
tion in these derivatives, and this analysis was 
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again verified by pilot comments. Doublet maneu- 
vers performed at 100 kn with artificial tail ice 
were accompanied by elevator force lightening and 
a strong "burble" feedback in the controls at high 
negative pitch rates. This characteristic was 
assumed due to the high tail downwash angle induced 
in the pitch over portion of the doublet maneuver, 
causing a significant amount of separated flow to 
occur on the horizontal tail plane. This occurred 
at a low peak load factor of approximately 0 g, 
which as mentioned earlier in this report, was 
erroneously input because the cockpit g meter read 
incorrectly. 

Conclusions 


Based on the analysis and results herein, two 
fundamental conclusions were drawn regarding the 
use of a Modified Stepwise Regression algorithm 
in estimating the effects of icing on aircraft 
stability and control characteristics: (1) The 

effects of the artificial ice shape attached to 
the tail were measurable as changes in the longi- 
tudinal stability and control derivatives for all 
flight conditions tested, and (2) that a multi- 
plicative relationship was identified between the 
ensemble standard error and the estimated standard 
error from the Modified Stepwise Regression rou- 
tine. The numerical value representing the ratio 
of the ensemble standard error, to the estimated 
standard error from the program when assigned to 
each control or stability derivative, are indica- 
tive of the level of confidence one should have 
in the estimated value of that parameter. 

The results of this flight test program indi- 
cate that MSR is a viable technique for determin- 
ing the effects of icing on the longitudinal 
stability and control characteristics of an air- 
craft. Application of this flight test and analy- 
sis technique should prove very useful for 
validating stability and control predictions from 
wind tunnels or computational methods. 

Appendix 
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along with the kinematic relations 

6 = q cos 4> - r sin <|> (A7) 

<j> = p + (q sin 4> + r cos 4>) tan e (A8) 


where u, v, and w are the velocity components; 

4> and 0 are the roll and pitch angles; g is 
the accelerationdue to gravity; p is the air den- 
sity; m is the airplane mass; lx, Iy, Iz and 
Ixz and the moments and product of inertia; S, b, 
and c are the wing area, span and mean aerody- 
namic chord; and a dot over a symbol represents 
the derivative with respect to time. The total 
airspeed V is given by 

V = (u 2 + v 2 + w 2 )l ^ 2 


Here, Eqs. (Al), (A3), (A5), and (A7) 
describe the longitudinal motion while Eqs. (A2), 
(A4) , (A6) , and (A8) describe the lateral motion 
of the airplane. The angles of attack and side- 
si ip are given by 

a = tan _1 (w/u) and |3 = sin _1 (v/V) 

The aerodynamic forces and moments are repre- 
sented by the coefficients C a (a = X, Y, Z, 1, m, 
or n). It is postulated that these coefficients 
can be written as a Taylor's series polynomial 
expansion about an equilibrium trim condition as 

i 

C = 0 a + V 0 a .Ax. (A9) 

a a,o Z-f a, 3 3 

U1 

where 0, A and 0, . are the unknown param- 
a,o a, 3 

eter. For a = X,Z, or m, the independent 
variables x^ are 


Postulated Models 

The six-degree-of-freedom equations of motion 
for the airplane are: 


a, ^ y , and 6e 

and for a = Y,l, or n, they are 


u = -qw + rv 
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Furthermore , 

Ax j “ x j (t) - 
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TABLE I. - PHYSICAL CHARACTERISTICS 
OF THE DEHAVILLAND DH-6 TWIN OTTER 


Characteristic 

Low 

High 

Geometric: 



Wing area, m 2 

39.02 

Wing span, m 

19.81 

Aspect ratio 


10 

MAC, m 

1.981 

Mass, kg 

4150 

4600 

Inertias, kg m 2 : 



Ixx 

21 279 

21 787 

{yy 

A ZZ 

30 000 

31 027 

44 986 

48 639 

Ixz 

1432 

1498 




TABLE II. - CHARACTERISTICS OF FLIGHT TEST INSTRUMENTATION 


[21 channel digital recording capability of the following 
parameters such that following units, ranges and 
resolutions are provided.] 


Parameter 

Range 

Resol- 

ution 

Longitudinal acceleration, g units 

*1 

0.002 

Laterial acceleration, g units 

±1 

.002 

Normal acceleration, g units 

+3,-1 

.004 

Pitch rate, deg/sec 

*15 

.029 

Roll rate, deg/sec 

±50 

.098 

Yaw rate, deg/sec 

±30 

.059 

Pitch angle, deg 

±45 

.088 

Roll angle, deg 

±90 

.176 

Alpha, deg 

+15,-10 

.024 

Left aileron position, deg 

+19,-16 

.034 

Right aileron position, deg 

+19,-16 

.034 

Elevator position, deg 

±20 

.039 

Rudder position, deg 

±18 

.035 

Beta, deg 

Battery A voltage, V 

±15 

.029 

Airspeed, m/sec 

0 to 85 

.083 

Air temperature, °C 
Voltage system reference 

-18 to 38 

.055 

Altitude, m 
Reference voltage 1 
Reference voltage 2 

-150 to 3000 

3.08 


TABLE III. - MEAN VALUES OF STABILITY AND CONTROL 
DERIVATIVES CALCULATED FROM 45 REPEAT MANEUVERS 
USING MODIFIED STEPWISE REGRESSION. 

[Flaps, 0°; 120 KIAS; baseline (no ice).] 


Parameter 

Mean value, 
6 

Standard errors 3 

Ratio, 

SE<e>/s<e> 

se<0> 

s(6) 

Cz a 

-5.66 

0.0493 

0.0169 

2.9 

Cz q 

-19.97 

1 .047 

.386 

2.7 

Czse 

-.608 

.0281 

.0122 

2.3 

Cm a 

-1 .31 

.0147 

.0191 

.8 

Cnriq 

-34.2 

.6444 

.4314 

1.5 

Cm^e 

-1 .74 

.0248 

.0131 

1.9 


a SE(9): Ensemble standard error; s(9): average program 

estimated standard error. 
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(b) FULL SPAN VIEW OF THE ARTIFICIAL ICE SHAPE. 

FIGURE 2. - ARTIFICIAL MODERATE GLAZE ICE ATTACHED TO THE TAIL OF THE ICING RESEARCH 
AIRCRAFT. 

ORIGINAL PAGE J5 

OF POOR QUALITY 


ORIGINAL PAGE IS 

OF POOR QUALITY 



FIGURE 3. - FLIGHT TEST NOSE BOOM SHOWING a AND P VANE INSTALLATION. 
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AZ, g UNITS q. DEG/SEC 0, DEG d. DEG ELEVATOR, DEG 







FIGURE 4. - CONTROL SURFACE INPUT AND AIRCRAFT RE- 
SPONSE FOR A TYPICAL "2-1" ELEVATOR DOUBLET. 
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INITIAL VELOCITY, kn 


120 


FIGURE 6. - PITCHING MOMENT COEFFICIENT DERIVATIVE 
ESTIMATES FOR 0° FLAP SHOWING THE EFFECTS OF TAIL 
ICING USING THE MSR TECHNIQUE. 20 ERROR BARS ARE 
FROM THE MSR PROGRAM. 






INITIAL VELOCITY, kn 


FIGURE 7. - VERTICAL FORCE COEFFICIENT DERIVATIVE ESTI- 
MATES FOR 10° FLAP SHOWING THE EFFECTS OF TAIL ICING 
USING THE MSR TECHNIQUE. 20 ERROR BARS ARE FROM THE 
MSR PROGRAM. 
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INITIAL VELOCITY. 


FIGURE 9. - 95 PERCENT PREDICTION INTERVALS ON ELEVATOR 
EFFECTIVENESS VALUES FOR THE ARTIFICIALLY ICED VERSUS 
UNICED AIRCRAFT TAIL. WING FLAPS AT 0°. 
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